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Abstract 

We examine an extension of the ideas of quantum cosmology and, in partic- 
ular, the proposal of Hartle and Hawking for the boundary conditions of the 
Universe, to models which incorporate Yang-Mills fields. Inhomogeneous per- 
turbations about a homogeneous, isotropic minisuperspace background model 
are considered, by expanding the Yang-Mills fields in harmonics of the spatial 
directions which are taken to be three-spheres. The expansions are made ex- 
plicit for SO(N) gauge fields thereby obtaining formulae compatible with the 
formalism conventionally used in quantum cosmology. We apply these results 
to the gauge group £0(3) and derive the Lagrangian and the semi-classical 
wave function for this special case. 



1. Introduction 



Recently we have witnessed a great activity in the interplay between high energy physics 
and cosmology The question of the initial conditions which in the early universe gave 
origin to the universe that we observe is clearly one of the main questions asking for 
clarification. Obviously, progress in such a direction needs an understanding of quantum 
gravity which is still lacking. On the other hand quantum cosmology provides a useful 
setting to discuss initial conditions and quantum gravity ideas. The most interesting 
results in this area were obtained from the so-called minisuperspace models where the 
infinite number of modes of the gravitational and matter fields was reduced by freezing out 
all the modes except of a finite number which are assumed to dominate . The quantization 
then proceeds through the application of the ADM [l[ formalism, thus obtaining the 
Wheeler-DeWitt equation as the Hamiltonian constraint and interpreting the conjugate 
momenta as operators, plus momentum constraints. It is clear that in order to go beyond 
the minisuperspace analysis one has to include the infinite number of the field modes, 
treated as frozen previously. This is usually done by treating these modes as perturbations 
around the minisuperspace degrees of freedom. 

In this context, Hartle and Hawking 0, |3| worked out the no-boundary proposal. An 
alternative suggestion is the so-called tunneling boundary condition 0]. Excellent reviews 
already exist on the subject ||. However, in most of the minisuperspace models and their 
extensions including perturbations which exist so far, the matter source is taken to be 
a scalar field f3[. Also fermion fields || and electromagnetism [|7| have been considered. 
Though this list contains fundamental fields, the absence of non-Abelian gauge fields 
is striking, especially if one considers that the evolution of the very early universe is 
dominated by gravitational and non-Abelian gauge fields. Recently a minisuperspace 
model including gauge fields was constructed || and the wave function for this model was 



determined || . A key ingredient of this model was the theory of symmetric fields [|10|, [IT 



which has been used before in the spontaneous compactification and dimensional reduction 
of higher dimensional models, to obtain SO (4)-symmetric ansatze for the fields in a Rx S 3 
topology. Let us remark here that it has been shown, at least in the case of a scalar 
field coupled to gravity, that if one considers perturbations of a minisuperspace model 
with three-spheres as spatial sections, any wave function, the object that we are mainly 
interested in, has to be SO (A) invariant [fl^| . 

The main subject of this paper is precisely to consider the above model beyond mini- 
superspace, including all the modes of the fields. For that purpose, we use the tool 
of harmonic analysis of symmetric fields on homogeneous spaces [11] , 13 , including to 
lowest non-trivial order the inhomogeneous modes of the gauge fields, while treating the 
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homogeneous degrees exactly. The method of harmonic expansion around a symmetric 
background configuration allows us to develop an analysis for arbitrary gauge group. 
Furthermore, and more importantly, it allows an analysis of the full realistic system of 
gauge fields coupled to fermions and scalars in a unified way. Therefore, using the method 
of symmetric fields one can consider the study of the coupled matter sources as well as 
supergravity theories coupled to Yang-Mills (YM), something which we consider as a 
continuation of our work. 

The organization of our paper is the following. In section 2 we review the construction 
of the symmetric model. In section 3 we apply the tools of the harmonic analysis on 
homogeneous spaces and symmetric fields on our particular system, i.e. a symmetric 
gauge field on R x S 3 with gauge group G = SO(N). Then we specialize on the case 
G = 50(3). In section 4 we present the perturbed model and in section 5 we consider 
the equations of the tensor, vector and scalar modes and their contribution to the wave 
function. We examine a particular case, namely the full system of perturbations, when 
the classical background field is fixed to its minimum value. Furthermore we discuss 
the interactions between gauge and gravitational perturbations, which might be quite 
interesting for the discussion of cosmological models and anisotropies in the cosmological 
background radiation. Finally in section 7 we present our conclusions. We also provide 
two appendices with details about the theory of symmetric fields and harmonic analysis 
on S 3 . 

2. General Formalism and the Symmetric Model 

Let us start by reviewing the relevant parts of the work already existing on the subject; 
this will allow us to introduce the problem and fix the notation. We will deal with closed 
universes of the form M = R x B, where B is a compact three-dimensional hypersurface 
which eventually will be the three-sphere; time t takes values in R. For the description of 
such a universe, a 3 + 1 formulation of the four dimensional Einstein theory is appropriate. 
Such a description, corresponding to an embedding of a three surface in the four dimen- 
sional spacetime is given by the ADM formalism [|(//, v = 0, . . . , 3, a, (3 = 1, . . . , 3) : 

^ = (D 
= -(N 2 - N a N a )dt 2 + 2N a dx a dt + h a pdx a dx? . 

The quantities N and N a are called lapse and shift functions, respectively, while h a p is 
a three-dimensional metric corresponding to B. The action of the Einstein- Yang-Mills 
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(EYM) system that we are going to consider has generally the form : 



S — S gr + Sym > (2) 



S, 



m 



GR 



2 r 

P 



167T 



/ d 4 x (-g)? (R-2X) + 2 / d 3 xh^K 
Jm JdM 



(3) 



SVm = A / Tr(F^F^) , (4) 

where R is the scalar curvature corresponding to g^ u , h = det (h a p) and if is the trace of 
the extrinsic curvature K a p of B = dM. It may be shown that eq. (|3]) may be brought 
to the form 

2 

Sgr = ^t f d 3 xdtN [ G a ^ s K af3 K jS + Vh( 3 R- 2A) ] , (5) 

1D7T J 

where 3 R is the curvature of B and G a ^ & is the DeWitt metric given by 

Q*fh6 = ly/h( h on/ h P6 + h a8 h fh _ ^ K <& ) . (6) 

The corresponding 3 + 1 - decomposition of eq. (f|) is 



Sym = —7 / d 3 xdt\fhNtY 
8e 2 J 



-2^F 0a F p - Ah^V^F^Fp, 



(7) 



One may proceed as usual to determine the conjugate momenta and the Hamiltonian. 
It turns out that the Hamiltonian is a sum of constraints : 

H = f d 3 x^h(N H + N a H a + tr(A H YM )) • (8) 

JdM 

The lapse and shift functions play the role of Lagrange multipliers, enforcing the con- 
straints which correspond to Einstein's equations, namely the momentum constraint H a 
and the Hamiltonian constraint H . Gauge invariance is guaranteed by the constraint 
H Y m with the time component A of the gauge field as the corresponding Lagrange mul- 
tiplier. 

Now we are ready to consider the quantized model, which will be described by a wave 
functional ty[h a p(x), A(x)] of the metric h a p(x) and the gauge field A(x). The constraints 
H , H a and Hym are expressed by demanding that their operator versions annihilate 
the wave functional ^\h a p(x),A{x)\. They ensure the invariance of the wave functional 
under reparametrizations and gauge transformations. The solution of the Wheeler-DeWitt 
equation H ^ = subject to the other constraints lies at the heart of the subject of 
quantum cosmology. 
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We proceed by reducing the infinite dimensional superspace (h a p(x) , A(x)) to a fi- 
nite dimensional set, the so called minisuperspace. Some of the very common defining 
conditions of minisuperspace are : 

1. The lapse is taken to be homogeneous, i.e. N = N(t). 

2. The shift N a is set to zero. 

3. The three-metric h a p is constrained to be homogeneous and isotropic. 

In quantum cosmology one is usually interested in closed universes. This, together with 
the above requirements fixes the spatial sections of M to be S 3 . Another way to single 
out this case is to impose 50(4) = SU(2) x SU(2) invariance on the metric of B = dM 
which leads to[| 

fi'oo = -a 2 N(t) 2 , g ab = h ab = a 2 a(t) 2 S ab , g 0a = . (9) 

An additional factor a 2 = 87r/m 2 P has been introduced for convenience. Using this metric 
we obtain 



Sgrq = — 3V53 / dt 



V -Na + H 2 Na 3 



N J ' < 10 > 

where H 2 = 87rA/3mp and V53 = 2tt 2 is the volume of the three sphere with radius one. 

For the matter sector, the following requirement is imposed, taking advantage of the 
gauge freedom : the change of the Yang-Mills field for motion on B at fixed t is allowed 
to be a gauge transformation rather than zero. In other words, we are imposing SO (4) = 
577(2) x SU{2) invariance on the gauge field modulo gauge transformation JET]. For the 
gauge group G = SO(N), according to the rules given in appendix A, the gauge potential 
is 

A(t) = A (t)dt + A a (t)cu a , (11) 

with 

Mt) = l^Wu, (12) 
A a (t) = [1 + X o(t)]T a + X \t)T aI , (13) 

and the Maurer-Cartan forms uj a on the coset S 3 = SU(2) satisfying du a + e a bc uj b Au c = 
0. The matrices Tab are representation matrices of SO(N) in correspondence with the 
decomposition SO (N) D SO (3) x SO(N - 3) 

N(N-l) , > , (JV-3)(JV-4), 

" = (3,1) + (I," ^ ") + (3,iV-3), 

I I (14) 

T a Tij T a i 



1 We will use Latin letters from the beginning of the alphabet to denote representation indices taking 
values in the orthonormal bundle. 
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which satisfy the following non-trivial commutation relations 



[Ta,T b ] 

[T a , T bK ] 
[Tij, T KL ] 
[Tij, T aK ] 



^-abcXcK i 

fijKT a i — SirTcJ ■ 



(15) 



The functions Xo(t), Xi(t) an d A km (t) are arbitrary. We should also note that the isomor- 
phic image of the isotropy group SO (3) in the gauge group SO(N) is taken to be the group 
generated by the three antisymmetric matrices Tj. We see in this example how the SO (4) 
symmetry requirement for the gauge fields has reduced an infinite number of degrees of 
freedom down to a fairly small number of functions. Moreover, for 5*0(4) symmetric fields 
the Lagrangian does not depend on the spatial coordinates, so the integration over the 
three sphere will just yield the volume V53 = 2tt 2 of S 3 as an overall factor. Taking this 
into account the action (|7|), after substituting 

F 0a = XoT a + 2x I T aI - 2 X iA IJ T aJ , (16) 
Fab = (1 - Xo - X 2 )tacbT c - 2xoe acb x I T cI (17) 



for the field strength becomes 

3 



5- 



YMo 



v s ,- 2 jdt 



a 

N 



,dxo, 
1 dt ' 



+ 2-V( Xo ,x) 
a 



where 



VtX = j t X + A%, 



V( X o,x) 



1 



(x^ + x 2 -i) 



;is) 



(19) 
(20) 



Given the action So = Sgr + Syai , we may find the conjugate momenta and the Hamil- 
tonian. It turns out that (dropping the factor 31^3 / e 2 for simplicity) 



Tin 



a . 
— a, 

N ' 



7T 



Xo 



X 



a ^ 

-P t Y. 

N A 



The Hamiltonian reads 



H = 



IN 
2~a 



2 2 1 tt2 4,2,2 
-7T a - a + H a + 7T X0 + 7T^ 



giving rise to the Wheeler-DeWitt equation 



1 
2 



d 







da \ da 



-a z + H 2 a 4 - 



d 2 



2 



dxl dx 2 



2V(Xo,x) 



2V( X0 ,X) 



4>(a, Xo,x) 



(21) 
(22) 

Eip(a,Xo,x) ■ 
(23) 
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The parameter p enters to take care of the operator ordering problem 0. For the case 
under investigation the gauge constraints H YM are related to the Lagrange multipliers A 
(cf. eq. (H)). They boil down to the classical condition 

Xi^xj ~ Xj^xi = °' ( 24 ) 

or, in the quantum language 

d d 

{Xia Xjt— ma, Xo, X) = 0. (25) 

OXJ OXI 

Based on the fact that the gravity and the gauge parts of the problem decouple, the 
classical behavior has been determined for a general SO(N) gauge group. It exhibits 
interesting features, for example wormhole solutions arise upon consideration of a Eu- 
clideanized version. Moreover, in the special case N = 3, the wave function has been 
determined in the semi-classical approximation by Bertolami and Mourao ||. Writing 
the separation ansatz 

y(a, X o)=J2 C n(a)Un(xo) , (26) 

n 

where C n (a) is the gravitational wave-function obtained by Hartle and Hawking [|J, 
U n (xo) is given by 

U( X o) = Aexp(-^(3( Xo + l) 2 -(Xo + l) 3 )) , for X o < -1 , 
U( X o) = 5exp(-^(3( X o + l) 2 -(Xo + l) 3 )) 

+ Bexp(-^(3(xo-l) 2 + (Xo-l) 3 )) , for | Xo | < 1 , 

U( X o) = Aexp(-^(3( X o-l) 2 + (xo-l) 3 )) , for X o > 1 ■ 

Here A = (1 + exp(— 4/3)) 5 and B is a normalization factor. The Hartle-Hawking 
boundary proposal has been imposed to arrive at these solutions by requiring that Xo( 7 7 — 
— oo) = =Fl. Furthermore Xo{v — 0) = Xo is the value at the given three-surface and rj is 
the conformal time. 



3. Harmonic Analysis 

A very interesting issue is to see what happens beyond the minisuperspace approximation. 
This has been done by Halliwell and Hawking for the case of a scalar field coupled to 
gravity. They consider perturbations around the minisuperspace model, considering the 
following forms for the various fields involved 

h a p = a 2 {Q a p + e a p), <P{y,t) = <P(t) + 5<P(y,t), N(y,t) = N (t) + 5N(y,t), (28) 
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where Q a p is the metric of the three sphere, allowing also N a (y, t) to take (small) non-zero 
values 0. The background quantities a, 0, N are treated to all orders, and the perturba- 
tions up to second order. The method to handle the perturbations is the use of harmonic 
expansion. In this paper we would like to examine the problem of perturbations around 
minisuperspace for the case of the SO (4) symmetric gauge fields referred to above. We 
therefore split the gauge field 

A = A (0) + A (29) 

into a symmetric part as it was presented for the group G = SO(N) in the last 
section, and a perturbation A. On the other hand, following appendix A, Aq (which 
behaves as a scalar as far as the SO (4) rotations are concerned) and A a can be written 
in the form 



A (t,y) = £ D^(L(y))a ^(t), (30) 

m,pq 

A a (t,y) = A B a + J2 D%\L(y))xa$f>(t). (31) 



rn.pq 



We note that the factors appearing in equations (|A.6|) and ( |A.7|) have been absorbed 
in the coefficients ao pq(t) and Xa ffi{t)- In these equations the field A^ can be identified 
with the background A B plus the ^/-independent parts of the expansions corresponding 
to the trivial representation (m = 0) of SO (4). In order to specify the harmonic analysis 
of the gauge field we must first solve the constraints given in appendix A for the 
5*0(4) representation matrices. 

We decompose the adjoint of the gauge group G under R = 50(3) ~ SU{2) according 
to the embedding : 

G D RxH , 

(32) 

adjG = (3,l) + (l,adj H) + J2 k (R k ,H k ) . 

Let us next introduce the generators T a , T L and T rkhk of G in correspondence with the 
decomposition (P^f). We write their nontrivial commutation relations as 



[T a ,Tb\ — e a b c T c , (33) 



[T a ,T rhhk ] = D^l{T a )T, hk . (34) 



The coefficients a j^'(t) and \a ^(t) take values in the Lie algebra of G and may 
be expanded in the basis of the matrices T g ,Tl and T Tk h k defined above. The relevant 
2 Here and in the following y denotes the coordinates of the three sphere. 
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expansions read : 



°o£>(*) = E«StoT c + E«S(^+ E ^lf h f k \t)T rkhk , (35) 

k,r k ,h k 



Xa%Ht) = ExSc(*) T c + ExS|l(*) T x+ E X^£ k \t) T nh k - (36) 

k,r k ,h k 



In (^) the indices r^, /i^ run over the representations Rk, Hk respectively. 

The above expansion has been made with an eye to the solution of the constraints 
( |A.8| ) and ( |A.9| ) in appendix A. Indeed, the constraints can be expressed in terms of 



the components of a ^'(t) and Xa P q{t) as follows (summation over repeated indices is 
implied) : 

D^KT c )a^(t) - e cba a^(t) = 0, 

Dff(T c )a™(t) = 0, (37) 

n( m ) (T \ / -,,( m )(- R fc'- ff fc) (+\ _ n( R k) (rp \ (m)(Rk,Hk) ( + \ _ rv 
^ps l J cJ a pq\r k h k \ Z ) LJ r' k r k \ 1 c) a sq\r' k h k \ l ) ~ U ' 



(3f 



* ) (T c )xS |J (t) + £ca6 xS, J (t) = o, 

U ps ( i cJ X a ,pq\r k h k + £ ™b X b ,sq\r k h k \ t ) 

T^{Rk)(rp\ {rn){R k ,H k ),,-, _ n 

U r' k h k \. 1 ^) X a ,sq\r' k h k \ l > - U. 

In this form it is easy to invoke Schur's lemma and find the solution of the constraints. 
According to the discussion in appendix B, we denote the representation (m) by (j i; j R \ J); 
in this notation the nontrivial solution is written as 



^lJrII) _ r (j L ,j R \l) 

_0i,ifl|0) _ (j L ,j R \0) 



(39) 



AjL,jR\J)(R k ,H k ) _ r AjL,jn\J)(R k ,H k ) -r /■ • \ — . 7~) 



X^sq\f ] = <la,lb\Js> X q jL ' jRlJ \J = 0,l,2, 

A3l,3r\1) _ r ('An") 

AjL,3R\J){Rk,H k ) _ p „ I T„ v. AjL,3R\J)(R k ,H k ) 

Xa,sq\r k h k ~ ^ 1Q > ^WftK 5 ^ Xg|^ 

In the last equation J takes the values jfc — l,jk,jk + 1) where jfc is the angular mo- 
mentum describing the representation R&. We conclude that the possible values of J 
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which determine the degree of harmonics which appear in the harmonic expansion are 
{0, 1,2, j k - 1, 3k, 3k + 1}- In particular: 

J max = m&x(2,j max + 1) , (41) 

where j max is the maximum angular momentum corresponding to the representations Rk 
which appears in the decomposition (^) of G under 50(3). The harmonic expansions 
read : 

+ Efc^o E q , L Dt M0 \mh%' jm W L (42) 

+ E {Rk , Hk) E q , ULM ^(^))4 Jalw ' Hk) m khk , 



A a (t,y) = A B a 

+ Ej=o I i 1 2E (i£liH )3jE 6jPl? ^^ |J) (^(y)) < i«, i6|Jp > x? ijH|J) (W 
+ E(i? fc ,H ft ) Ej=ij fc -i,ii fe ,ii fc+ i E(,- £ j«)3j Ep l(? ,r fc ,h fc D p q L,m]{J \ L (y)) 

^L^R^iRkHk) (4-\rp 



<la,R k r k \Jp> X % L J R]JmkHk) (t)T ri 



(43) 

The expression D^ q L ' jRl0) (L(y)) on the second line of (g) just means that the first index of 
this harmonic takes the value corresponding to the identity representation in the product 
(3l,3r)- On the other hand the notation J2(j L ,j R )3j means that we sum over all the 
tensor product representations (j'l, J_r) of 5*0(4) w SU{2) L ® SU(2) R , which contain the 
representation J upon restriction to SU(2)ai ag , as described in appendix B. 

If the gauge group G is SO(N), several simplifications occur. The decomposition of 
the adjoint of SO(N) under R = 50(3) can be written in a more concrete form and it 
is given in eq. (0). The commutation relations of the generators, heavily used in the 
calculations, corresponding to the decomposition ( |T4"D can be found in eqs. (|T5|). 

Maybe the most important simplification brought about is that the only representation 
appearing in the decomposition (T4) is the triplet, corresponding to angular momentum 
jmax — 1- This makes things much easier, since according to the statements made above 
only scalar- vector- and tensor-harmonics appear in the expansion. In this case we are 
also in a position to give explicit expressions for the relevant Clebsch-Gordan coefficients. 
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Changing from the SU(2) standard basis to the "tensor" basis they read 



< la, 1&|00 > -> 8 ab 

<la,16|lp> -> e abp (44) 

< la, 16|2 Pxp 2 > -> K^apx^bpj + <WA P2 ) _ |^a^p lP2 • 

All indices run from 1 to 3 and we have represented the J = 2 components by the double 
index pip 2 , as in Appendix B. 

We are going to discuss in the following section the case N = 3. Most of the physical 
features of the model are present in this case; on the other hand, calculations are much 
easier, because of some additional simplifications. Also the formulae to be presented, 
although far from being simple, are much more economical. The simplifications coming 
from choosing N = 3 are summarized below. 

One may see that the parts of the expansion having to do with adj SO(N — 3) and 
(Rk,Hk) (cf. eq. (|32|) ) no longer exist. Thus the harmonic expansion becomes : 

Mt,v) = E^h)^!^^* 11 ^^))"?^ 11 ^*)^, (45) 

A a (t,y) = A^ + Ej=o^E {jL j R ^jE b , P , q D^ J \L(y)) 



< la,lb\Jp > xi jL ' jRlJ) (t)Ti 



(46) 
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We separate the background field and the y-independent part in the last equation and 
write explicitly the sum over J : 

A a {t,y) = {l+ Xo {t))T a 

+ Eo i ^ )3 oE M ^ jB|0) (^(y))ta^^ |0) (t)T fc 

+ E 0i ^) 3 iE^^^ |1 H^(y))^ P x^^ |1) Wr 6 (47) 

+ E(j L ,j R )D2 J2b,q,p!,p2 [<j + <5&pi$ap 2 ) 

- ^a^ PlP2 ]^( ( ^ 2) (^(y))X^^ |2) (i)T6, 

where we have plugged in the explicit expressions (041) for the Clebsch-Gordan coefficients. 
The Clebsch-Gordan coefficient \{S api 8b P2 + Sb Pl S ap2 ) — \8 a b8 PlP2 in fact coincides with the 
operator F '(pif^ltfift) > which projects out the part of the harmonic D^'J^\L(y)) which 
is symmetric and traceless in the indices p\ and p 2 . Thus, this Clebsch-Gordan coefficient 
just projects out the harmonics G ( ^\L(y)) } S^\L(y)) and P^ q (L(y)) } discussed in 
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Appendix B, where it is also explained that the harmonics D^ q L ' : ' R ^ 1 '(L(y)) will give rise to 
S^\L{y)) and P$(L(y)). Similarly, the harmonics D^ Mo) {L{y)) become Q q n \L(y)). 
Then, introducing some numerical factors for later convenience, the final form of the 
harmonic expansion in this case reads 



A {t,y) 



EnE b , q [V2a^\t)S^ + \L(y)) + s/2af-\t)S§-\L{y)) 



+ T^tKt)P^{L{y))\T b 



(4? 



A a (t,y) = 4 0) 



with 



+ 



7 ^(t)5 ab Q^(L(y))T b 



EnE 6 , 9 j= 2 P q n+) (t)S% g + \L(y)) + ^r\t)SX~\L{y)) 



in-). 



+ E n E b , q ^ n) (t)P^(L(y))]e acb T b (49) 
+ E n n, q [^ +) (t)GX>(L(y)) +nt-\t)G^{L{y))]T b 

+ £nE& !<z V6[€jp\t)P$ q (L(y))]T b , 

4 0) = (l + Xo(t))T a . (50) 



4. The Perturbed Model 

Let us now consider inhomogeneous perturbations around our minisuperspace. We use 
the exact action in terms of the background quantities but expand only to second order in 
the perturbations. We assume that the metric is of the form given in eq. (||) except that 
it is multiplied with the normalization factor a 2 . The gravitational part of our model is 
treated in exactly the same way as in ref. 0. Therefore let us recall that 



h ab = a e a (5 ab + e ab ) 



(51) 



where e a = a and (in simplified notation) 

v/6 



S ab a n Q n + V6 b n P: b + y/2 c n S n ab + 2 d n G n ab . 



(52) 



Furthermore for the lapse and shift functions we have : 

1 



N 



Nn 



v/6 



9nQ n 



(53) 
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N„. 



(54) 



Let us remark here that, of course, these expansions are exactly the same as the ones we 
get, if we apply the symmetric field constraints on the functions e a b, N and N a . Now the 
gravitational part of the action (|5[) takes the form B 



Sgr = Sgr + V53 ^2 / dtLQ R 



(55) 



where 

^GR 



~-e a N 



le 3a 

+ 2iVo 



3V - \)< + ^ trA^l - 2 (n 2 - 4) c 2 



n- 



1 



2 2 
-{n 2 + l)^ + -(n 2 - 4)a n b n + -{n 2 - 4)g n b n 

2 , 2 l ^ 

+ -)g n an 



1 



3(n 2 - 1) 

n 2 - 4 • 



/c 2 + (n 2 - 4)j 2 



al 4 



< "Or, 



n- 



b z n + (n 2 - A)ci + d 2 r: 



2,CL n (l n -{ 



n 



+ 8(n 2 - A)c n c n + 8d n d n + 2g 



n 



b 2 n + 6(n 2 - 4)c 2 + 6d 2 n + 3a n g n - g 2 n 



-3e 3a Nn H 2 



+e~ a 

1 , 



-a. 



4 



2 n z — 1 

, / 2 2n 2 -4- 2 \ . 2 
K I _ 3 a « ~~ 3 n 2 _ + 3 a #« I ~ 2 ( n ~~ 4)c n j n 

_ n ZiA h 2 <~ 2 ^ 2 
n 2 -l " 



(n - 4)c n - d n + -a n #„ 



(56) 

For the YM part of the action with gauge group G = 5*0(3) we write the field strength 



as 



F a 
F n ,h 



(o) 



F^ + V^A b -V^A a +[A a ,A 



(o) 



(57) 



with the symmetric field given in eq. ( |17|) and a covariant derivative with respect to 
the symmetric field : 

Vi 0) = V a + [4°),-] , VS 0) =«9 . (58) 



Inserting the symmetric field ( |i~3"| ) and the harmonic expansions (f48|) , (f49|) together with 
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the above expressions in eq. (0) we obtain^ 



with 



and 



Sym — Sym + -Jj- X! / dtL YM , (59) 



L" M = -^£J + — C n 2 + rr^^ + — Q (60) 
yAi e° 1 N 2 N e a 3 iV 



^2,^ 4 2 3 2 3 . 2 i \ 2 



^n 2 - 4 
n 2 — 1 



+Xo(Xo - 1) ( a nln + -&n£n + 2(n 2 - 4)c n z/ n + 2d n p n - g n ^ ri 



2 / 3 2 2 1 1 2 1 2 n 2 -4, 

-(n 2 - 4)c n p n - -#„cr n + n(ra 2 - 4)c n z/„ + 2nd n p n 

( 1 ~ 1 / 2 

+Xo I -7nOn + -npnpn - np n p n - -n(n - A)v n u n 

- 4)7* - \(n 2 - 4K " g^f^ " ~(- 2 - 4)nV 2 
1 2 1 2 n 2 — 4 

+ l )^n ~ gfa ~ 4 )£n - + 3 ^2 _ 1 hn ° n + ~~ 4 ) C «P« 

1 2 1 

+ 3^^ - ™(™ 2 - 4)c n z/ n - 2nd n p n - -(n 2 - 4)7„£„ + -n{n 2 - 4)p n v n , 

(61) 



• 2 / "^2 X 77/ ^72 ~^ / 2 ,i \ 2 ~^ ?2 2 ^ ^ 

+X ( ~o a n + n3 A + o( n - 4 ) C n + o rf n + 7#n ~ 7°^, 



2 n 2 - 1 " 2 v ' " 2 " 4"" 4 
/ 1 1 2 n 2 — 4 n 2 — 4 

+ X0 -dnPn + -SnTn + 'KPn ~ 4— -b n £ n + 2{r? - 4)c n Q! n 

6 2 3n - 1 n z — 1 



-2(n 2 - A)c n i) n - 2d n p n - ~g n f3 n ~ \g n in + 2a n p n + - 1 3 n a n 

6 2 Sn z — I 

I n 2 — 3 i 2n 2 - 4 

+ (n 2 - 2)a 2 - (n 2 - A)a n u n + — -f3 2 n + ~/? n7n - g^Ti^ 

+ 7^ 2 + \pn + \ 2 1 ; °n + \pn + \{n 2 - A)i) n 2 + 2^- — ~ na n p n , 
4 2 on - 1 2 2 rr — 1 

(62) 

A = (x^-l) 2 f 6(n2 1 _ 1) fc n + 2 ^) > (63) 



3 Here we give the Minkowskian action. For the Euclidean version the signs of and £3 have to be 
reversed. 
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A = XoXo (~3 ^J3 i fcnCr » ~~ 2 -^) +Xo(Xo - !) (~3 n 2 _ i kn @ n ~ 4 in«n) 
/ 1 2 n 2 — 4 ~ \ 

+ X0 ( -g^n7n - 3 2 _ ~ ^ 2 ~~ 4 )i" Z/ n + J ( 64 ) 

In these (and the following) formulae a sum over the even and the odd part of the relevant 
coefficients is implied. The tilde over a coefficient signals an off-diagonal term connect- 
ing even and odd perturbations with each other p(^ ven /° dd ) — p(pdd/even) f or a g ener j c 
coefficient p n . 

One may ask about the remnants of gauge symmetry in this Lagrangian. To answer 
this question we start with an ordinary infinitesimal gauge transformation SA^ = V M £7 — 
[U, Ap], U = u a T a and demand that the symmetric field A^ remains fixed. The gauge 
transformation will therefore exclusively act on the expansion A and we find 5A^ = V^U 
in the leading order. Expanding the parameter Ui 

u a = -J=t; n P; + V2w n S: (65) 



and using the expression ( |i~3|) for the symmetric field as well as eqs. ( fl8|) and (49) for A 
this transformation property can be translated to the coefficients : 



(66) 



5u n = w n , 5a n = w n , Sp n = 2xoW n + nw n , Sfx n = 0. 
Our procedure is analogous to the definition of coordinate transformations acting on the 



perturbations of a certain background metric as given in ref. WM. This analogy can be 



pushed one step further by introducing gauge invariant variables 

T„ = 7n + 2£ n , B n = (3 n + 3j n , S n = a n + 3xo7n , , . 

(67) 

A n = a n - i> n R n = p n - 2x Q v n ~ nv n . 

As in the case of coordinate transformations in pure gravity the tensor mode /i is invariant, 
since the expansion (B5I) does not contain any tensor degree of freedom. It might be 



This mixing of the odd and even parts may give rise to questions about parity conservation when 
the perturbations are included. This does not create any problem though: If we start from the parity 
transformation properties of the metric and the gauge fields, we may derive the ones of No(t), a(i), Xo(t) 
and the perturbation coefficients. Using these properties, it turns out that the second-order Lagrangian 
is parity invariant. 
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interesting to note f\ that the above variables represent a linearized version of the invariant 
quantities tr(B c B d ) introduced by Liischer |L5| (B c = e ac bF ab is the magnetic field.). With 
these variables one gets, for instance: 

V2 



Rr. 



dfi(y)tr(B c B d )S: d . 



(6* 



8(n 2_ 4)(x 2_ 1) 

The important observation is now that our Lagrangian (|60|) is invariant under the 
transformations fl66|) provided the field Xo fulfills the background equation of motion to 
be derived from the action (|l"8|). Consequently, in such a case the Lagrangian can be 
written in terms of the gauge invariant variables only. For a special solution xo we will 
demonstrate this explicitly in section 5. 

One can proceed from this point by defining the conjugate momenta and obtaining 
the Hamiltonian in the usual manner. Here we want to indicate the structure of the 
Hamiltonian 



H = A^o + ^^+^n^ + Et^ . (69) 



fin) 



in) 



/(") 



which we find in accordance with eq. (^) . The linear YM gauge constraints and Tip 
responsible for the gauge invariance of the wave function read 



(n) 



T-/( n ) = 

nf = 

with the conjugate momenta 



2Xovr p „ - gP„7r x + ir Un + nn Pn 



i 



9(n 2 -l) an7l X 3"7> 



1, 



(70) 



7T V 



71% 



7T„ 



Tin 



Tin 



7T„ 



3 ^'n + (quadratic terms) , 



2e 2 N 

Vg*. 



In + hPn + Xo{an - 9n) 



e 2 N 



[Pn ~ 2xo"n - na n ] + -^{xl - l)jn 



3e 2 Vo(n 2 -l) \Pn X0A1] + 3e 27v („ 2 -l) (Xo ty^r. 

[fi n - 2xod n ] , 

n 



e 2 N 

V s3 e a (n 2 -4) r. 
e 2 N \rn 

4V" s 3e Q (n 2 -4) 
e 2 Af (n 2 -l) 



(71) 



Their quantized version can be used to generate our gauge transformations (p6|) per- 
forming the commutator 5A a = [A a , w n 7i^ + v n H^} with the harmonic expansion (^9|) 
5 We thank C. Lee for pointing out this relation. 
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for A a . The pure gravitational parts of eq. (|69|) and the momenta for the gravitational 
perturbations can be found in ref. ||. We will not give the other YM parts explicitly 
since we are mainly interested in the semi-classical wave function in a region of super- 
space where the minisuperspace variables are classical but the perturbations are quantum 
mechanical. Therefore, it is enough in our approximation to solve first the Euclidean 
background equations by setting the perturbations equal to zero and then solve the equa- 
tions of motion for the perturbations subject to the Hartle-Hawking boundary condition 
by neglecting their back reactions on the background. In this way, the wave function 
is given by \& = Cexp(— Ie) where C is a semi-classical prefactor and Ie the extremal 
Euclidean action. I E is given by Ie = (l/2)p n 7r Pn where p n is a generic perturbation. 



5. The Wave Function 

Let us follow the above program for the particular case of the EYM system which we are 
considering here. The background equations of motion for the minisuperspace degrees of 
freedom a, xo are given by 

e< (£) + ^ + W - S W + ^ + %V)=0, (72) 



^i(^rXo)+V = o, (73) 



dt\N J 2 N 2 u 2 u 2e 2 \2iV u e 

e a d (e a 
No"dl ViVo 
while for N we obtain the constraint 

c I i „-2a tj2\ 1 I e -2 , iV o 



2 \Nj +e -^-^nA^^V' - <74) 

Here V denotes the derivative of V = 1/2(1 — Xo) 2 (cf. eq. (p0[)) with respect to xo- 
The perturbations of the gravitational and YM fields can be separated into three types, 
namely tensor (d n , fi n ), vector (c n , u n , j n , a n , p n ) and scalar (a n , b n , a n , j n , f n , (3 n , g n , k n ) 
modes. In the equations of motion different types of perturbations do not couple. This 
observation allows to consider separately each type of perturbation and calculate their 
corresponding contribution to the wave function. 

For the tensor perturbations d n and \x n we obtain the following equations of motion 



^2X0 - Vj dn + V'fin + 2nUfi r 



e 2a 



-2j^X0^n 



v 



(75) 

~N~ Q dl {~No~^ n ) + ^ 2 + = 2 Jp^° dn ~ 2n ^°^« + 2nU - Vd n , (76) 



e a d ( e a \ , . e 



2d 
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with the abbreviation U = Xo ~ 1- 

Let us now consider the case of the vector perturbations. To simplify the calculation 
one can make a choice of gauge. For the gravitational field we choose to proceed in the 
transverse-traceless gauge which corresponds for the vector modes in c n = 0. For the YM 
field inspection of the transformations ( |66D shows that we are free to take p n = 0. In this 
gauge the equations of motion and the constraints of the remaining variables (V n , j n , a n ) 
are given by : 



e a d 

iVo It \N 



e a (n 2 -4) + 



1 2 



V 



e 2 e a 



2(n 2 + 2U)a n + 4nxo«r 



-XoKVn ~ TTfJnXO , (77) 
iV 

(2V'a n + X o{n 2 -A)u n + 2nUa n ) , (78) 

2(V'j n + nUj n ) + e a (n 2 - 4)z> n . (79) 



In the case of the scalar modes the Lagrange multipliers are k n , g n , f3 n leaving five 
genuine variables. The transverse traceless gauge for the gravitational degrees of freedom 
dictates for the scalar modes a n = b n = and we will choose o n = for the YM modes. 
The remaining modes (j n , £ n , /3 n , g n , k n ) fulfill : 



0(£<* + 5«) + l< n ' + B + ^ + l< B '- 4 >«- 



£s(^-5«) + i(" , - 1 ^ + 26 ') 



,2a 



j-^Xogn - 2V'g n 



2e c 



3NJ 



JO) 

ID 



Pn 



9 



1 



n 2 - 1 + 2U 



A{n 2 - 4)4 - {n 2 - l) 7n + {n 2 - l) X ^n + e~ a V'k r , 



n ~ 1 ■ r 



n 9 • 9n 

e ae 



1 3N 2 3 

2Xo(37n + Pn) + ^^'Tn - ^Xl9r, 

n 2 — 4 11 

X07n + 2— -Xo^n + ~— -V'fir, 



+ 3ae a g n 
2 1 



Vr — V 



2n 2 - 1 



+ 



1 



1 



ae a n 2 — 1 



52) 



-/Cr,. . 



e a n 2 - 1 

To calculate the wave function we need the value of the Euclidean action for the classical 
solution which in our gauge is given by : 

^y-e 2a (d n d' n + 4c/d 2 ) + ^ (// n // n - 2x' d n fx n ) , 

(n 2 - 4)(z> n - ct„> n , 



d 



2e 2 N 
Vs°e° 
e 2 N n 



3) 



1 1 n 2 — 4 • 1 

7(7n + - fin ~ Xo9n)ln + 2— -(£„ - -/3 n )£ r , 



rr 



6 
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In the above expressions the prime denotes the derivative with respect to the conformal 
time r] with fig = -iga. 

We proceed to solve the above system of equations in the large n approximation, 
i. e. we only consider the terms with the highest occurring power of n. After eliminating 
the auxiliary modes the equations decouple in this approximation and can be written in 



terms of the gauge invariant variables introduced in eq. (67). They take the simple form 
P" = n 2 P n for P n = d n , fi n , R n ,T n . Imposing the Hartle-Hawking boundary condition 
P n — > for i] — > — oo singles out the solutions P n ~ exp(n?7). With the eqs. fl83D this 
results in the wave function 





= C exp 


Vs3 e 2a ndl 
2 


Vs 3 2 




= C' exp 


Vs3 d2 _ 








= C'exp 


a 2 nT n 







We see that the perturbations start out in their ground state. In contrast to the gravita- 
tional mode d n the frequencies corresponding to the YM perturbations do not depend on 
the scale factor e a . This behavior is a consequence of conformal invariance of YM theory. 



A case of particular interest is when the classical field Xo obtains its minimum value 
Xo = ±1- Obviously this represents a solution of the background equation of motion (f73|) 
compatible with the Hartle-Hawking boundary condition. As already discussed, the La- 
grangian in such a case can be explicitly expressed in terms of the gauge invariant vari- 
ables. An additional simplification occurs because Xo — ±1 corresponds to a symmetric 
field which is pure gauge as can be seen from the eqs. (|16|), (|17j). This causes the 
decoupling of gravitational and gauge perturbations in the quadratic terms we are con- 
sidering here. The Lagrangian Ly M = T C n + v C n + s C n then takes the form : 



V C n 



N, 







R z n - 4xoR n A n - 2nR n A n + 2n 2 A 2 n + A Xo nA n A n 



2N 

4e a 



2N 
N 



Qe a 



n 2 Rl 

n 2 - 4 
n 2 — 1 

V 



2nx R n R, 
1 



55) 



f 2 + 

n 1 



Q2 

3(n 2 - 1) n 



2{n 2 - 4) 
3(n 2 - 1) 



Br,T r 



3(n 



^——XoBnSn + ^B, 



4)r 2 + s 2 

) n 1 n 
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Starting from these expressions it is easy to calculate the equations of motion and elim- 
inate the Lagrange multipliers A n and B n . Solutions can be found without any further 
approximation and we obtain the explicitly gauge invariant wave function 



with 



exp 
exp 

exp 



2e 2 

4e 2 

_Vs 1 
6e 2 



{nR 2 n - 2 X0 RnRn 



(r n ,s n )M^ 



M (n) 



( (n 2 + 2)(n 2 -A) - Xo (n 2 - 4) ^ 



n (n* 



n 



J 



*6) 



17) 



1 \ -Xo(n 2 - 4) 

As a consistency check one can apply the linear gauge constraints H^, and one finds 
that they annihilate and s ^ ( - n \ 



Finally we want to discuss some features of the general solutions to the non-trivial 
systems of differential equations which govern the classical evolution of the perturbations. 
The solutions considered above have the drawback that they describe situations, where 
the gravitational and gauge sectors are completely decoupled as for the symmetric con- 
figuration. The separation was due to the fact that the YM Lagrangian is scale invariant 
and the only symmetric gravitational degree of freedom is that of the scale factor. In the 
first case the n 2 -approximation is the first term of an approximation in rapidly changing 
fluctuations and the decoupling appears as a consequence of the diagonality of the kinetic 
terms 0. The splitting in the case where Xo is fixed to the minimum of the potential 
is again a property of the YM Lagrangian, together with the restriction to terms up to 
second order in the non-symmetric fields: From the general from of the YM action each 
term containing gravitational and gauge perturbations has to be supplemented by a com- 
ponent of the background field strength. It is therefore multiplied with a time derivative 
of Xo or a term proportional to Xo ~ 1 which both vanish for xo — ±1- This effect is an 
intrinsic property of the symmetric YM potential and its knowledge in the formalism used 
in this paper is crucial as opposed to the approach, where the YM degrees of freedom are 
treated from the beginning as scalar fields. Obviously the above reasonings do not apply 



6 In fact this approximation reduces to a large momentum one for large enough n due to local Lorentz 
invariance. 
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to expressions of higher order in the perturbations and therefore there is no reason that 
the decoupling survives even in these special cases. 

To investigate the interactions between the gauge and gravitational perturbations 
let us consider the next-to-leading-order equations in a large n approximation of the 
coupled system describing the tensor modes, which are relevant for the discussion of 
density perturbations. Invoking the ansatz 

dt =1/2(4— )± 4°^)) , /4 =i/2(fi^ en) ±^ dd) ) 



e n VD ± = e n VM ± 



n i 



we find the following set of equations at order n: 

-^-L>±-tanh(?7-^) J D± + ^cosh(7 7 -?7 a ) 2 / ± M± = 0, (89) 
ar] e z 

±M±+g ± M± + h ± D± = 0, (90) 

where dt] = iNsa~ 1 dt. We recall at this point that the exponential e nr> in equation 
is the leading order behaviour found previously, while the prefactors and 



represent the next-to-leading-order corrections. In eqs. the conformal mode a has 
been replaced by the classical solution 

gT 1 ^) = Hcosh(i] - r] a ) , (91) 

where r\ a is defined by the value of a on the given three-surface, a(0) = a . Inspection 
of equations (|89|) , (|90D reveals that the functions and do not really depend on 
n; from now on we simply denote them by D±, M±. Inspection of the same equations 
shows that, unlike the leading order exponential, the prefactors depend on the value of the 
field Xo(v)- This field will be replaced by the classical solution given in ref. 0, expressed 
through the functions /, g and h defined below: 

/± = ±(^(77) - 1) + x'oiv) , = ±Xo(v) , h± = ±(xl(v) - 1) - XM ■ (92) 

There are four different solutions, determined by their boundary values at rj = oo and 
7] = : 



(93) 





Xo(v) 


= coth(?7 — 


Vo) , 


Xo 


< -1 


i7 


Xo(v) 


= tanh(?7 — 


Vo) , 


\Xo\ 


< 1, 




Xoiv) 


= — tanh(r] 


~Vo) , 


\Xo\ 


< 1, 


IV 


Xo(v) 


= — coth(r/ 


-Vo) , 


Xo 


> 1 . 
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On the other hand, 770 expresses the boundary condition on Xoiji) at 77 = : Xo(0) = Xo- 
For instance, for the first of the above solutions, one gets: coth(— r] ) = Xo] this defines 
rjo for this solution. We should remark that Xo{v) wn l remain in the region from which it 
started off at 77 = 0. The solutions of eqs. ( jS9l) corresponding to the four possible choices 
for Xo{v) are 

I : 

D+ = A + cosh (77 - 7] a ) , 

M + = A + [cosh(?7o - 7] a ) sinh(r] - r/o) _1 + 2 sinh(r7 - T) a ) cosh(/7 - 770)] 

+B + sinh(r/ - 770) , ( 94 ) 

D- = -^-A- [cosh(r/o - r) a ) cosh(r/ - r) a ) sinh(?7 - 7? )~ 2 

+2 sinh(r/ — r) a ) cosh(r/ — r) a ) coth(r/ — 770)] + B_ cosh(^ — r) a ) , 

M_ = A_ sinh(?7 — 77o) _1 , 

II : 

D + = A + cosh(?7 - rja) , 

M + = A + [2 cosh (770 - T} a ) sinh(r7 - 770) 

- sinh(77o - 77 a ) cosh(77 - 770) _1 ] + B + cosh(77 - r] ) , (95) 
D- = -fr^4- [2 cosh(77o - 77 a ) cosh(77 - r] a ) tanh(77 - 77 ) 

- sinh(77 — r] a ) cosh(77 — rj a ) cosh(77 — r] )~ 2 } + £L cosh(77 — r] a ) , 
M_ = v4_ cosh (77 — r/o) _1 , 

III : 

D + = ^-A + [2 cosh(77 - 77o) cosh(77 - r] a ) tanh(77 - 770) 

- sinh(?7 - r) a ) cosh(7? - r] a ) cosh(?7 - r] o y 2 ] + B + cosh(?7 - 7? a ) , 

M+ = cosh(7? - r/o)" 1 , (96) 

D- = A_ cosh(7] — 77 a ) , 

M_ = — A_ [2 cosh (770 — 77 a ) sinh(77 — 77 ) 

- sinh(?7o - r] a ) cosh(77 - r/ ) _1 ] + 5- cosh(?7 - 77 ) , 
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IV : 

D+ = t§-A + [cosh (770- 77,,) cosh^-r/Jsinh^-rft)- 2 

+2 sinh(r] - ?? a ) cosh(?7 - r] a ) coth(r? - r/ )] + B + cosh(r/ - ?y a ) , 

M+ = A+sinh^-^o) -1 , ( 97 ) 

.D_ = A_ cosh(?] — ?y a ) , 

M_ = — A_ [2 sinh(?7o — r] a ) cosh (77 — r] ) 

+ cosh(?/ - ?7a) sinh(?7 - rfa) -1 ] + 5- sinh(7? - 770) , 

where A± and -B-t are integration constants. Note that there are two types of relations 
between coefficients entering both the gravitational and gauge perturbations: in the first 
case the relative strength depends on the ratio H/e while it is universal in the second 
one. Insofar as the physical properties of these solutions are concerned, it is interesting 
to note that the interactions can cause quite strong correlations between the two kinds 
of fluctuations. For instance, if the ratio ^ is normalized to a fixed value for a given 
conformal time r], it is predicted for any other value of 77 and can differ substantially 
from the initial choice. This might give rise to interesting consistency conditions for 
cosmological models, if this ratio is constrained for physical reasons at two different times 
of the evolution. 

Starting from the general expression ([83]) for the classical action and eliminating the 
derivatives via the eqs. ( 8^ , 90) one finds the tensor wave function in the 0(n) approxi- 
mation : 



V n) = Cexp 

_Vs 1 
2e 2 



e 2a {n + 3a')d 2 n (9? 



+ Xo^nfin + 2(1 - xl)d n Hn ~ ^Xo d nf^r. 



(99) 



As expected it shows a mixing between the gravitational and the YM mode in the new 
subleading terms. 



7. Conclusions and Outlook 

In this work we have extended the minisuperspace approximation treated in previous 
works to include the full infinite-dimensional model of symmetric gauge fields coupled to 
gravity. A primary ansatz characterizing our approach is that we treated gauge fields, such 
that the change of the gauge potential due to the translation on S 3 can be compensated 
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by a suitable gauge transformation. 50( TV) -symmetric gauge fields have been considered, 
however the extension to unitary, symplectic or exceptional groups is within reach; in fact 
it should be made, if a realistic model is to be constructed within this framework. 

We restricted the gauge group to be just 5*0(3), to keep the formulae tractable. We 
feel that this restriction does not hide too many characteristics of the system. Most 
probably, a general SO(N) group would have richer topological structure, but we have 
not addressed this issue in this paper. Probably working with different class of gauge 
groups would give rise to essentially different phenomena. So, e. g., the groups SU(N) 
with a conventional embedding of 5*0(3) would impose the use of spinor harmonics. 

Only a part of the superspace has been investigated, namely the portion where the 
minisuperspace variables are classical, while the perturbations are quantum mechanical. 
To investigate other regions, one should try to solve the Wheeler-DeWitt equation in 
a better approximation. This would allow taking into account the back reaction of the 
perturbations on the background fields. An interesting feature of our calculation has been 
the existence of gauge invariant variables and the possibility to express the second order 
Lagrangian in terms of them, illustrated in the main text for the case xo — ±1- 

We will present the detailed predictions of the wave function found above in a future 
publication. For the time being we comment on some very basic characteristics of this 
solution. The wave function has been found in the regime of large n (the order of the 
corresponding harmonic). 

As a first step, we only kept the terms of order n 2 ; the wave function we found shows 
that the perturbations start off in their ground state. This means that approximating the 
inhomogeneous and anisotropic degrees of freedom just to second order is well justified 
and describes the system quite well in the relevant region of the superspace. We also 
point out that the vector perturbations are genuine in our case and not a gauge artifact, 
as is the case of the scalar fields. 

When we proceed with the next-to-leading order, i. e. we also keep terms of order n, 
we encounter new, quite interesting, phenomena. Let us pay particular attention to the 
tensor perturbation modes, which are e. g. relevant for the spectrum of relic gravitons in 
the universe. We stress that, in contrast to models with scalar fields, we have a second 
tensor mode in the theory stemming from the YM field. In the next to leading order this 
mode couples to the gravity tensor perturbation often referred to as the linear graviton. 
A similiar behaviour shows up in the corresponding tensor mode wave function. It is 
conceivable that these effects could have interesting cosmological consequences. 

One line of development of our work would be the inclusion of scalar and fermionic 
fields and consider the full system at work. 
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Moreover, one may calculate wormhole effects using our results. For example one 
can find the effective interactions between gauge fields at large distances caused by the 
presence of wormholes. This requires calculating the matrix element of the product of 
two gauge field operators between the flat space vacuum state and wormhole states. The 
latter one can be computed semiclassically starting from our equations of motion for the 
perturbations and finding solutions which fulfill the boundary conditions appropriate for 
wormhole states. 
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Appendices 

A Symmetric Fields 

Let us briefly summarize here how one can perform the harmonic analysis of symmetric 
fields on homogeneous spaces. For the general method of constructing symmetric fields 



we refer the reader to refs. [1 LJ where a detailed treatment is given. Let us just comment, 
that a field is called S symmetric on a space M x if a symmetry transformation by 
an element of the isometry group S corresponds to a gauge transformation. This implies 
that functions of the fields such as the energy-momentum tensor or the Lagrangian density 
are independent of the coordinates of S/R just because they are gauge invariant. The 
requirement that transformations of the fields under the action of the symmetry group 
of S/R are compensated by gauge transformations, leads to certain constraints on the 
fields. The solution of these constraints provides the theory on M. In particular, a gauge 
field A M on M D splits into A m on M and A a on S/R; A m behaves as a scalar under S- 
transformations and lies in the adjoint representation of the gauge group G. The surviving 
gauge symmetry H on M is that subgroup of G which commutes with R. In other words 
the gauge group H in one dimension is the centralizer of R in G, i.e. H = C G (R a ); Rq 
is the isomorphic image of R in G. The remaining components of the gauge field A a 
become vectors under the coset space transformations. The transformation properties of 
the fields A a under H can be found if we express the adjoint irreducible representation of 
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G in terms of Rq x H 



G D Rq x H , 

(A.l) 

adjG = (adj/U) + (l,adj#) + £,(i4^). 

and 5 under R 

S D R, 

(A.2) 

adjS 1 = adji? + J2i Si ■ 

Then for every pair R iy Si where R4 and Si are identical irreducible representations of R, 
there remains a multiplet on M transforming under the representation Hi of H. All other 
scalar fields vanish. 

In order to perform the harmonic expansion of the fields on M x S/R we have to find 
a suitable ground state configuration around which to expand. Since the symmetry of the 
vacuum M x S/R is the symmetry group of M times S, we demand that the background 
configuration Af, should possess these symmetries. Therefore has to vanish while A B 
has to be S-symmetric. Such a configuration is given by 

K = eU , (A.3) 

where a = 1, . . . , dimS 1 — dimi?, i = 1, . . . , dinrR, e z a are vielbeins of the coset space and Jj 
are the generators of the gauge group G spanning the algebra of R G . The background field 
A% is also S'-symmetric in the broader sense i.e. including a suitable gauge transformation. 
Indeed we obtain that 

5 A A B a = -(d a {Y A )Ji + [O^Jj, A B a ] , (A.4) 

where Q\ are the so called i?-compensators, and therefore the gauge transformation 
needed to compensate this S'-transformation is given by : 

W A = -Q\J t . (A.5) 

Let us now discuss the harmonic expansion of a field on a coset space. It is well known 
that the matrix elements of the inequivalent unitary irreducible representations of 
S, indexed by (m) , serve as a complete orthonormal basis into which any function defined 
on S can be expanded (this is known as the Peter- Weyl theorem in group theory ||16|| ). 
Since S is the isometry group of S/ R, we can still use, in the case of the coset space, the 
matrix elements of the unitary irreducible representations of S as an orthonormal basis 



for the harmonic analysis, although these now have to be properly restricted |T7[]. The 
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harmonic expansions of the fields A m , A a can be written as 



A m {x,y) 

A a (x,y) 



m' Pa V «X> 



X 



m pq 




u>1> 



A B a + A a 



(A.6) 
(A.7) 



where L(y) is a representative element of each inequivalence class, d m is the dimension of 
the representation m of S and dv is the dimension of the T> representation of R. As we will 
see, only those irreducible representations of S contribute to the harmonic expansion of a 
field on S/ R in the representation T> of R, which contain T> in their decomposition under R. 
Furthermore, all the relevant D^^s contribute to the harmonic expansion only with those 
of their rows which correspond to the particular /^-representation T>. In order to determine 
the relevant representations of S which contribute to the harmonic expansions ( |A.6| ), 
( |A.7| ), we should demand that under a general ^-transformation, which now includes 
the gauge transformation given in eq. ( |A.5| ), the coefficients in these expansions should 



transform according to some irreducible representation of S JT^]. We find 



(m) 
pq 



+ fiab4>b 



M _ r j. 

sq IVf 



M _ 

m pq 



[Ji: 



7 Hi 

*m S q J 



0. 
0. 



(A.8) 
(A.9) 



Eqs. ( |A.8[ ), ( |A.9[ ) specify which representations of S are present in the harmonic expansions 
of the corresponding fields given in eqs. (A.6), ( ] A . 7| ) . From the solution of eq. (|A.8 ) we 
find that the relevant representations of S for the harmonic expansion of A a are those 
which include in their restriction to R the representations of R appearing in the cross 
product of Ri and Si given in eqs. (|A.l ) and ( |A.2| ). Similarly, for the harmonic expansion 
of A m one should include only those representations of R which appear in eq. ( |A.1| ). All 
these 5* representations will contribute only with their rows which correspond to the R 
representation. Furthermore, the ^-symmetric fields are the ^/-independent ones, i.e. the 
first terms in the expansions ( |A.6| ), ( |A.7| ). 



B Harmonics 



Various versions of the harmonics on S 3 have appeared in the literature [B5|, [T^, 2D|. 
We use the ones basically proposed in [O, which bear great similarity to the harmonics 
used by Halliwell and Hawking, but their properties are much more transparent and easy 
to derive from group theoretical considerations. We refer the reader to [IH for a more 
complete account and present here some basic steps, needed to adapt these considerations 
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to our notation. The manifold is S 3 , or SU(2), which is considered as a coset space of 
SU{2) L ® SU{2) R over SU(2) diag , where SU{2) diag is defined as usual as the subgroup 
of SU(2)l ® SU(2)r with equal left and right rotations. The harmonics for a field on 
the above coset corresponding to an angular momentum J, say, will transform as the 
irreducible representations of SU(2) diag . Such representations may be sieved out of 
the tensor product representations of the group SU(2)l <8> SU(2)r which contain 

upon restriction to SU(2) d i ag . The tensor representations are characterized by the 
two angular momenta ji,jR corresponding to SU(2)l, SU{2)r respectively. For instance 
the scalar harmonics (corresponding to J = 0) are contained in 

(0,0), (§,§), - , (B.l) 

while the vector harmonics (J = 1) are contained in 

(o, l), (i, 0), (l, l), (i |), (|, |), ... , (B.2) 

As a matter of notation, we denote the harmonic characterized by J constructed out of 
the pair (jL,jii), by D^ q L '^ R][J ^{L{y)). The indices p and q label matrix elements and 
L(y) is the representative element for the coset, at the point described by y in some local 
coordinate system. The additional index £ is used to distinguish between the various 
£>( J )'s that may be contained in D^ m ' . In fact for the SU(2) group only one appears 
for fixed J, so we completely drop this index in the sequel. We also note that the index p 
is constrained to run only over the 2 J + 1 values corresponding to the relevant T>^ J \ As 
explained in Appendix A, any field on the coset space transforming according to T>^ can 
be expanded as follows : 

4 J) (L(y))= £ J2^' jR]J)D fc' jR]J) ( L (y)) ■ (B.3) 

The orthonormality relation reads : 

E J S3 MV) [Dt jR{J \L{yW [D^ J \L{y))\ = ^J+i^ + i) ^^^ ■ 

(B.4) 

As the reader will notice, we define in the sequel the transverse harmonics with a normal- 
ization equal to y // ( 2 ^+ 1 j( 2 J^+ 1 ) _ motivation was simply to make the right hand side 

of eq. (|B.4| ) equal Vs^ L ^5^ R ^8 qq i . It is useful to define the quantities 

n = j L + 3r + 1, d = j L - j R , (B.5) 
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in terms of which the Casimir operator of the (jl, jr) representation of 50(4) SU (2) L <g> 
SU{2) R equals 

C 2 (jL,jR)=n 2 + d 2 -l. (B.6) 

Two important results read : 

V a D^JH\J){ L (y)) = - \T^\ r D%^\ J )(L(y)) , 

V a V" Dfc>i«\ J \L(y)) = [C 2 (VW)-C 2 (j L J R )}D^x\ J \L(y)) (B.7) 

= [J(J + 1) - (n 2 + <? - 1)] Dfe^l^)(L(y)) . 

In the above the T^ L ' iR ) are the part of the reductive decomposition (T^> L ^ R \T^ L ' jR) ) of 
the generators of (jljJr), which corresponds to SU(2) diag ; C 2 (T>^) = J (J + 1) is the 
Casimir operator for the representation J defining the harmonic. These important results 
mean that any relation involving differentials on S 3 reduces to mere algebra. Two more 
identities are needed to find the properties of the harmonics. These are the commutation 
relations for covariant derivatives acting on vectors and tensors : 

(V a V c -V c V a )H = V d R d bac , 

(B.8) 

(V a V 6 - V h V a )H cd = H ed Rt ab + H ce R dab , 

Rabcd = hJ>ad — &achd- (B-9) 

In the following we give explicitly some properties of the scalar, vector and tensor har- 
monics, which are the only relevant ones for the case G = SO(N), whose special case 
N = 3 is considered in this work. 
1) Scalar harmonics 

These are scalar eigenfunctions of the Laplacian. It is convenient to define them with 
the following normalization : 

Q^ ] (L(y))=n 2 D^°\L(y)). (B.10) 

The first index on DQ q L ^ R ^\L{y)) takes the value corresponding to the identity represen- 
tation, so there is no point in keeping it in the notation for the harmonic Q„ n ^- We also 
use the index (n) defined above to characterize the scalar harmonics, rather than the pair 
{JliJr)- I n fact for the scalar harmonics, Jl = ]r = j, so the other combination, namely 
d, vanishes in this case. The quantity j runs over integer and half-odd-integer values. For 
the normalization introduced above, it is easy to show that 

J s3 dn(y) iQ q n) (L(y))Y [Q^'\ L {v))\ = v s ^'5 qq , . (b.h) 
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The eigenvalue of the Laplacian is found from the difference of the Casimir operators 
(eq. flBJ) : 

V a V a Q q n \L(y)) = - {n 2 - I) Q q n \L(y)) . (B.12) 

2) Vector harmonics 

The vector harmonics on a simply connected coset space without boundary can be 
uniquely decomposed into a sum of the "transverse vector harmonics", characterized by 
the vanishing of their covariant divergence and the "longitudinal vector harmonics" , given 
by the gradient of some (scalar) function. This corresponds to the fact, obvious from 
eq. (|B.2|) , that the quantity d may take three values, namely —1,0, and +1. Thus the 
vector harmonics split into categories, according to the value of d. The ones correspond- 
ing to d = —1 (d = +1) are called odd (even) transverse vector harmonics; we denote 
them by S^~\L(y)), S^ q + \L(y)) respectively, while for d = we get the longitudinal 
vector harmonics, P^j(L(y)). The index a runs from 1 to 3. We note that the terms 
"odd" ( "even" ) derive from the behavior of the corresponding harmonics under a parity 
transformation [jnj. Let us now be more specific and give the exact definitions of the 
various harmonics : 



Sj$t\Hv)) = \j^D^\L(y)). (B.13) 

On the other hand, since the scalar harmonics, previously defined, span the space of scalar 
functions on the coset space, the longitudinal vector harmonic may be chosen to be the 
gradient of a scalar harmonic : 



P$(HV)) = ^ZiV a Qf\L{y)). (B.14) 



Using the above definitions, we obtain the following properties of the vector harmonics : 

E a S s s d»(y) [S^\L(y))\* [S^\L(y))\ = VgsS^S* , 

EaIs»Mv) [S^iHvW [SSf/\Hy))] = 0, (B.15) 

v a sSf (%)) = o , 

V a P$(L(y)) = -Q q n KL(y)), (B.16) 

Sa, b V a S^\L(y)) = nS^\L(y)), 
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V«V* (L(j/)) = -(nZ-2)S^\L(y)), 
V„V" P^Lfo)) = -(n 2 - 3) P^Lfo)) . 



3) Tensor harmonics 

As in the previous case, also the tensor harmonics can be classified into several kinds, 
according to the value of d = ]l — jr. This quantity may now take the values ±2, ±1, 0. 

The harmonics with d = —2 (d = +2) are called odd (even) transverse tensor harmon- 
ics. To get them from the D^^ r]{2 \ we represent the index p (which may take 5 values) 
by a pair of indices aia 2 , each running from 1 to 3. Then the definition is 



G^ q (L(y)) = E PM qi q 2 )D^f(L(y)), (B.18) 



where the definition of the projection operator P(aia 2 \qiq2) reads : 

P(aia 2 |gi<?2) = [2(^101^202 + <W<W) - g<W^oia 2 ]- (B.19) 

This operator projects out the symmetric and traceless part of the expression it acts on. 
Thus the harmonics G^^ q (L(y)) are symmetric, traceless and transverse. We find the 
following additional properties of the transverse tensor harmonics : 



(B.20) 



(B.21) 



E ai , a J^(y)[G'itU^(y))]1<c:U L (^)] = °» 

^*C(£(y)) = ~(n 2 - 3)G^ k (L(y)) , 

^GSWy)) = 0, 

^CiWf)) = > 

^«V ttl G5,(L(l,)) = nG^Lfo)). 
The remaining types of tensor harmonics can be expressed in the form : 

E P(aia 2 |gig 2 )V gi D^ |1) (L(y)). (B.22) 

91,92 

Basically they are covariant derivatives of harmonics with smaller "magnetic quantum 
numbers" d, properly symmetrized. We recall that any vector harmonic is a linear com- 
bination of the transverse and longitudinal vector harmonics. Thus these remaining har- 
monics are linear combinations of : 
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1) The even and odd "longitudinal-transverse" harmonics, with \d\ = \]l ~ 3r\ = 1- 
Eq. flBgp becomes : 



£ (»)) = 2[ V - S - ±)(L(?/))+V -^ ±)(L(?/))] (B - 23) 

and 

2) the "longitudinal-longitudinal" tensor harmonics, with d = ji — jn = 0. Now 
eq. (B.22|) becomes 



P^UiHv)) = ^aAtim) + Va 2 P^ q (L(y))) - l -8 ai a 2 V a P^{m) , (B.24) 
respectively. From the definition of P^(L(y)) one can prove that 

PSZl\ q (Hv)) = ^r[V ai V a2 QP(L(y)) + l -5 aia2 Q^\L{y)\ (B.25) 

which is the definition of this quantity used by Halliwell and Hawking. 
In the following we gather some properties of these harmonics. 

Ea 1A /dMy)[^ ) | 9 (L(y))]*[5^( L (y))l = v^^^-s^s^ , 

Ea lia J^(y)[^S g (i(y))]lCl|V( L ^))] = °. (B.26) 



3(^-1) 

^S { :H(L(y)) = -(n a -4)^(L(y)) 

,(n±) , 



V ai V«^S 9 (L(y)) = 0, 

v B v.52K(^(y)) = -(^ 2 - 6)^:3,(^(2/)) , 

V a V a Pi;i 2k (L(y)) = -(n 2 - 7)P£l [q (L(y)) , (B.27) 

V ai V»P£l ]9 {Hv)) = §(n 2 -4)QW(L(y)). 
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